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Abstract 

The Johnson-Lindenstrauss lemma allows dimension reduction on real vectors with low distortion 
on their pairwise Euclidean distances. This result is often used in algorithms such as fc-means or k 
nearest neighbours since they only use Euclidean distances, and has sometimes been used in optimiza¬ 
tion algorithms involving the minimization of Euclidean distances. In this paper we introduce a first 
attempt at using this lemma in the context of feasibility problems in linear and integer programming, 
which cannot be expressed only in function of Euclidean distances. 


1 Introduction 


In machine learning theory there is a wonderful and deep result from functional analysis and geometric 
convexity called the Johnson-Lindenstrauss Lemma (JLL) [3]. Intuitively, this lemma employs a concen¬ 
tration of measure [3] argument to prove that a “cloud” of high-dimensional points can be projected in 
a much lower dimensional space whilst keeping Euclidean distances approximately the same. Although 
this result was previously exploited in purely Euclidean distance based algorithms such as fc-means [1] 
and k nearest neighbours [3] (among others), it has not often been applied to optimization problems. 
There are a few exceptions, namely high dimensional linear regression [5], where the application of the 
JLL is reasonably natural. 

In this paper we present some new results on the application of the JLL to establish the feasibility 
of Linear Programs (LP) and Integer Linear Programs (ILP). We consider problems with m constraints 
(where m is large), and reduce m by projection on a random subspace, while ensuring that, with high 
probability, the reformulated problem is feasible if and only if the original problem is feasible. 

The geometrical intuition underlining our idea stems from the cone interpretation of LP feasibility. 
Let P be a feasibility-only LP in standard form (i.e. all inequalities have been turned into equations by 
the introduction of m additional non-negative variables), written as Ax = h with a; > 0, where A is an 
m X n rational matrix, b S K"* is a rational vector, and a; is a vector of n decision variables (which might 
be continuous or integer). Then P can be interpreted as the following geometric decision problem: given 
a cone spanned by the columns of A, is b € cone(A) or not? In this setting, the role of the JLL is seen 
to be the following: if we project A and 6 in a smaller dimensional space, the JLL assures us that the 
“shape” of the projected cone and of the ray b are approximately the same, and hence that the answer 
to the problem will be the same most of the times. 

In Section 2, we formally define the problem. In Section 3, we recall the Johnson-Lindenstrauss 
Lemma and prove some results linked to its application to the ILP case. In Section 4, we derive results 
concerning LP feasibility when the cone generated by the matrix A is pointed. In Section 5, we generalize 
the previous results, proving that the distance between a point and a closed set should be approximately 
preserved. Finally in Section 6, we present some computational results. 

^Supported by a Microsft Research Ph.D. fellowship. 
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2 Linear programming and the cone membership problem 


It is well-known that any linear program can be reduced (via an easy bisection argument) to LP feasibility, 
defined as follows: 

Linear Feasibility Problem (LFP). Given b G K™ and A G Decide whether there 

exists X G K" such that Ax = b A x > 0. 

We assume that m and n are very large integer numbers. We also assume that A is full row-rank. In 
particular, we have m < n, since otherwise we can find x uniquely from Ax = 6 by taking the left inverse 
oi A. 

LFP problems can obviously be solved using the simplex method. Despite the fact that simplex 
methods are often very efficient in practice, there are instances for which the methods run in exponential 
time. On the other hand, polynomial time algorithms such as interior point methods are known to scale 
poorly, in practice, on several classes of instances. 

If oi,..., a„ are the column vectors of A, then the LFP is equivalent to finding a: > 0 such that b is 
a non-negative linear combination of ai,... ,a„. In other words, the LFP is equivalent to the following 
cone membership problem: 


Cone Membership (CM). Given &, oi,..., a„ G K™, decide whether b G conejai,..., a„}. 


This problem can be viewed as a special case of the restricted linear membership problem, which is 
defined as follows: 

Restricted Linear Membership (RLM). Given 6, ai,...,a„ G and X C R”, decide 

n 

whether b G linjc(oi,..., a„), i.e. whether 3A G X s.t. b = XiUi. 

i=l 

For example, when X = R", we have the cone membership problem; and when X = IP (or {0,1}") we 
have the integer (binary) cone membership problem. 

It is known from the JLL (see below for an exact statement) that there is a (linear) mapping T : 
R™ ^ R^, where k m, such that the pairwise distances between all vector pairs {ai,aj) undergo low 
distortion. In other words, the complete graph on {ai,...,a„} weighted with the pairwise Euclidean 
distances realized in R™ can also be approximately realized in R^. We are now stipulating that such 
a graph is a reasonable representation of the intuitive notion of “shape”. Under this hypothesis, it is 
reasonable to expect that the image of G = cone(ai,..., a„) under T has approximately the same shape 
as C. Thus, given an instance of CM, we expect to be able to “approximately solve” a much smaller 
(randomly projected) instance instead. Notice that since CM is a decision problem, “approximately” 
really refers to a randomized algorithm which is successful with high probability. 

Notationwise, every norm || • || is Euclidean unless otherwise specihed, and we shall denote by E'^ the 
complement of an event E. 


3 Random projections and RLM problems 


The JLL is stated as follows: 
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3.1 Theorem (Johnson-Lindenstrauss Lemma [3j) 

Given e € (0,1) and A = {ai,...,a„} be a set of n points in R"*. Then there exists a mapping 
T : R™ —>■ R^ such that 


(1 -e)||a* -a^ll < \\T{a^) -T{aj)\\ < (1 + e)||a, - aj|| 


for all 1 < *, J < n, in which k is 0{e ^ logn). 


( 1 ) 


Thus, all sets of n points can be projected to a subspace having dimension logarithmic in n (and, 
surprisingly, independent of m), such that no distance is distorted by more than 1 + 2e. The JLL is 
a consequence of a general property (see Lemma [3.21 below 1 of sub-Gaussian random mappings T{x) = 



P X where P is an appropriately chosen matrix. Some of the most popular are: 


• orthogonal projections on a random fc-dimensional linear subspace of R™; 

• random kxm matrices with each entry independently drawn from the standard normal distribution 

A^(0,1); 

• random kxm matrices with each entry independently taking values +1 and —1, each with proba¬ 
bility 

• random kxm matrices with entries independently taking values -1-1, 0, —1, respectively with 
probability 1, |, 1. 

3.2 Lemma (Random projection lemma) 

Let T : R™ —>■ R^ be one of the above random mappings. Then for all e € (0,1) and all vector x € R™, 
we have: 

Prob( (1 - e)||x|| < ||T(x)|| < (1 + ^)||x||) > 1 - (2) 

for some constant C > 0 (independent of m, k,e). 

Not only can this lemma prove the existence of a mapping satisfying conditions in the Johnson-Lindenstrauss 
lemma, but also it implies that the probability of finding such a mapping is very large. Indeed, from the 
random projection lemma, the probability that Eq. dl]) holds for a\\ i ^ j <m is at least 

1 _ 2 Q = i-m(m- (3) 

Therefore, if we want this probability to be large than, say 99.9%, then simply choose any k such that 
ioom(^m-i) ^ This means k can be chosen to be fc = which is 0(£“^(ln(m) -|- 

3.5)). 

We shall also need a squared version of the random projection lemma. 

3.3 Lemma (Random projection lemma, squared version) 

Let T : R™ —R^' be one of the random mappings in Lemma Td.21 Then for all e € (0,1) and all vector 
X € R™, we have: 


Prob( (1 - e)||a;f < \\T(x)\\^ < (1 + £)||a:f ) > 1 - (4) 

for some constant C > 0 (independent of m, k,e). 

Another direct consequence of the random projection lemma is the concentration around zero of the 
involved random linear projection kernel. 




3 RANDOM PROJECTIONS AND RLM PROBLEMS 


4 


3.4 Corollary 

Let T : R™ —>■ be one of the random mappings as in Lemma 13.21 and 0 ^ x G . Then we have 

Prob(T(x) ^ 0) > 1 - 2e-^^ (5) 

for some constant C > 0 (independent of n, k). 

Proof. For any e G (0,1) , we define the following events: 

A = {T(a:)^0} 

B = {(l-e)|N|<||T(x)||<(l + e)||x||}. 

By Lemma 13.21 it follows that Prob(S) > 1 — 2e“°^ ^ for some constant C > 0 independent of m, fc,e:. 
On the other hand, .4*^ fl S = 0, since otherwise, there is a mapping Ti such that Ti{x) = 0 and 
(1 — e)||a;|| < ||Ti(a;)||, which altogether imply that x = 0 (a contradiction). Therefore, B C A, and we 
have Prob(.4) > Prob(S) > 1 — This holds for all 0 < e < 1, so Prob(4l) > 1 — □ □ 


3.5 Lemma 

Let T : K™ —>■ be one of the random mappings as in Lemmaand 6, oi,..., a„ G R"*. Then for any 

given vector x G R", we have: 


(i) Ifb= Y, XiQi then T[b) = Y XiT{ai); 

i=l i=l 


(a) If b ^ Yi=i Prob 


Tib)^Ytl^^na^) 


> 1 - 2e 


-Ck. 


(in) Ifb^ Sr=i all y G X C R", where |X| is finite, then 


Prob 


Vy GX r(6) 


i=l 


> 1 -2|X|e 


-Ck. 


for some constant C > 0 (independent of n, k). 

Proof. Point jH) follows by linearity of T, and ([ii| by applying Cor. 13.41 to Ax — b. For (lull) , we have 


Prob 


VyGXT(6)^^y,r(a,) 


Prob 






= 1 - Prob 






2=1 
[by (E 


v^x 


i=l 


> 1 - ^ = 1 - 2|X|f 


-Ck 


y&x 


as claimed. 


□ 


□ 


This lemma can be used to solve the RLM problem when the cardinality of the restricted set X is 
bounded by a polynomial in n. In particular, if |X| < n'^, where d is small w.r.t. n, then 

Prob[r(6) ^ Linjf {r(ai),...,r(a„)}] > 1 - 2n‘'e-‘''=. (6) 

Then by taking any k such that k > ^ ln(|) + ^ Inn, we obtain a probability of success of at least 1 — 5. 
We give an example to illustrate that such a bound for |2f| is natural in many different settings. 
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3.6 Example _ 

If X = {x G {0, 1}"| X]r=i some d, where 0 < for all 1 < i < n, then \X\ < n‘^, where 

n n n _ 

a*; then ^ ^ which implies ^ d.- 

i—1 i—1 i—1 

Therefore |X| < (g) + (")+... < n^, as claimed. □ 

Lemma [S3] also gives us an indication as to why estimating the probability that 

T{b) ^ cone{r(ai),..., r(a„)} 

is not straightforward. This event can be written as an intersection of infinitely many sub-events 

n 

{r(6)^5^y,T(a0} 

i=l 

where y € M" ; even if each of these occurs with high probability, their intersection might still be small. As 
these events are dependent, however, we still hope to find to find a useful estimation for this probability. 


d = max 1-^1. To see this, let a = min 

l<i<n “• ~ l<i<n 


4 Projections of separating hyperplanes 


In this section we show that if a hyperplane separates a point x from a closed and convex set C, then its 
image under a random projection T is also likely to separate T{x) from T{C). The separating hyperplane 
theorem applied to cones can be stated as follows. 

4.1 Theorem (Separating hyperplane theorem) 

Given b ^ conejai,..., a„} where 6, ai,...,a„ £ R™. Then there is c G R™ such that c^b < 0 and 
c^’^Ui > 0 for alii = 1,... ,n. 

For simplicity, we will first work with pointed cone. Recall that a cone C is called pointed if and 
only if C n —C = {0}. The associated separating hyperplane theorem is obtained by replacing all > 
inequalities by strict ones. Without loss of generality, we can assume that ||c|| = 1. From this theorem, 
it immediately follows that there is a positive £o such that c^b < —Eq and c'^Oi > Eq for all 1 < f < n. 

4.2 Proposition 

Given b,ai,... ,an G R™ of norms 1 such that b ^ cone{ai,..., a„}, £ > 0, and c G R™ with ||c|| = 1 be 
such that (Fb < —e and c^Ui > e for all 1 < i < n. Let T : R'” —>■ R^ be one of the random mappings as 
in Lemma \3.3[ then 

Prob[T(6) ^ cone{T(ai),..., T(a„)}] > 1 — 4(n -|- -G)k 

for some constant C (independent of m, n, k, e). 

Proof Let A be the event that both (1 — £)||c — x\\^ < \\T{c — x)\\'^ < (1 + £)||c — a;|p and (1 — 
£)||c -I- cclp < ||T(c -I- a;)|p < (1 -f £)||c -I- a;|p hold for all x G {6, oi,..., a„}. By Lemma 13.31 we have 
Prob(A) > 1 — 4(n -|- For any random mapping T such that A occurs, we have 

= i(||T(c+&)f-||T(c-&f) 

< J(l|c+&f-||c-bf) + |(||c + &f+ ||c-6|n 
= c^b -I- £ < 0 


{T{c),T{b)) 
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and, for all i = 1,..., n, we can similarly derive c^Ui — e > 0 from {T(c), T{ai)). Therefore, by Thm. 14.11 
r(6) ^ cone{r(ai),... ,r(a„)}. □ □ 


From this proposition, it follows that the larger e will provide us a better probability. The largest e 
can be found by solving the following optimization problem. 

Separating Coefficient Problem (SCP). 

Given b ^ cone {ai,..., a„}, find e = max {e| e > 0, c^b < —e, > e}. 

c,e 

Note that e can be extremely small when the cone C generated by oi,..., a„ is almost non-pointed, i.e. 
the convex hull of oi,..., a„ contains a point close to 0. Indeed, for any convex combination x = 
with Ai = 1 of Oi’s, we have: 


n n 

ll^^ll = ll^^ll l|c|| > XiS = e. 

i=l i=l 


Therefore, e < min{||x|| | x € conv{ai,... ,a„}}. 


5 Projection of minimum distance 


In this section we show that if the distance between a point x and a closed set is positive, it remains 
positive with high probability after applying a random projection. First, we consider the following 
problem. 

Convex Hull Membership (CHM). 

Given &, oi,..., a„ G R™, decide whether b G convjai,..., a„}. 

5.1 Proposition 

Given ai,...,a„ G R™, let C = convjai,..., a„}, b G R'” such that b ^ C, d = min||6 — a;|| and 

x£C 

D = max ||6 — ai\\. Let T : R™ —R*^ he a random mapping as in Lemma Id. 21 Then 

l<i<n 

Prob[T(6) ^ T{C)] > 1 - (7) 

for some constant C (independent of m, n, k,d, D) and e < 


Proof. Let Sg be the event that both (1 — e:)||x —i/lp < \\T{x — y)\\'^ < (l + £)||a: — ?/|p and (1 — e)||x + 
< 11 ^( 2 ; + ?/)|P < (1 + e)||2; + 2/|P hold for all a;, ?/ G {0,6 — oi,..., 6 — a„}. Assume Sg occurs. Then 
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n 

for all real > 0 with ^ = 1, we have: 

i=l 


||r(b) - XiT{ai)f = II XiT{b - m 


^A?||r( 6 -a 0 ll '+2 ^ NXj{T{b-ai),Tib-aj)) 

i=l l<i<j<n 


- (1 + e) fli - aj 


= f^A?||r(6-a0ll' + 5 E A,Ah||r(&-a, + b-a,)||^-||r(ai-a,)||A 

i—1 l<i<,j<n ^ A 

n y 

^ (1-£)E^?II^““*II^ +2 ^ XiXj - e)\\b - ai + b - ajW^ - (1 

i=l l<i<j<n 

n.n . 

= ll^-E^*®‘ll^ “^( + 2 ^ AiAj(||b - fli+ fo-O jII^ + ||ai - ajll^) j 

i—1 A l<i<ij<n A 

n ^ n \ 

= Wb-'^XiCif - e(^X'^\\b-atf + E + 11^ “ ) ■ 

i—1 A l<i<_j<7i ^ 

From the definitions of d and D^ we have: 

n , n /nx.2 

||r(6)-E^ir(ai)||^ >d2-eD^( E^? + 2 E AiA,) = - eZ)^ f E ) = 

2 ^ 

= d?—eD'^ > 0 due to the choice of e < Now, since ||r(&)—X]> 0 for all choices of A, it fol¬ 
iar 

lows that T{h) ^ conv{r(ai),... ,T(a„)}. In summary, if occurs, then T[b) ^ conv{r(ai),... ,T(a„)}. 
Thus, by Lemma l3.3l and the union bound. 


Prob(r(fe) i T{C)) > Prob(5'e) > 1 - 2(n + 2 ( ^ ) )e“ 


= 1 - 2ne~ 


for some constant C > 0. 


□ 


□ 


In order to deal with the CM problem, we consider the A-norm of x S cone{ai,... ,a„} as ||xm = 

n n 

min { X! I — OAx = ^ Xiai}. For each x G conejai,..., an}, we say that A G M!} yields a minimal A- 

i=l 

n 

representation of x if and only if ^ Ai = HxH^. We define /ta = inax{||x||A | x S conejai,..., a„}A||x|| < 

i—1 

1}; then, for all x G conejai,... ,a„}, ||x|| < ||x||a < ^A||a;||- In particular > I- Note that p-A serves 
as a measure of worst-case distortion when we move from Euclidean to || • ||a norm. 

5.2 Theorem 

Given 5, ai,...,a„ G M'" of norms 1 such that b ^ C = conejai,..., a„}, let d = min||6 — x|| and 

x£C 

T : R"* ^ be one of the random mappings in Lemma 1,3.31 Then: 

Prob( T{b) i cone{r(ai),..., T(a„)} ) > I - 2n(n -P i)e-CC'-'^')fe (8) 

for some constant C (independent of m, n, k,d), in which e = j^2^p\\iiA+i ' 

Proof. For any 0 < e < I, let Sg be the event that both (1 —e)\\x — y\\'^ < ||r(x —y)|p < (l-P£)||x — y|p 
and (1 — e)||x -P j/p < ||r(x -P ?/)|p < (I -P e)||x -P j/|p hold for all x,y G {b, ai,..., a„}. By Lemma [33J 
we have 


Prob(S'e) > I - 4^"" 2 j = 1 - 2n(n -P 
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for some constant C (independent of m,n,k,d). We will prove that if Se occurs, then we have T(b) (ji 
cone{T(ai),..., T(a„)}. Assume that occurs. Consider an arbitrary x € conejai,..., a„} and let 

n 

^ Xiai be the minimal ^-representation of x. Then we have: 

2=1 


||r(fo)-r(a;)||^ = ||r(6)-^A.r(a, 


= ||T(fo)f+ ^A?||r(aOf-2^A,(r(6),r(aO)+2 ^ XAi{T{ai),T{aj)) 

i=l i=l l<i<.j<n 

= l|r(6)f+^A?||r(a,)f+^^(||T(6-aOII 




> (1 - £)ii6f + (1 - £)E + E - “^11' - (1 + ")ii^+“^11') 


+ E 


^^((1 - e)\\ai + Ujf - (1 + e)\\ai - ajf), 


because of the assumption that S'g occurs. Since ||&|| = ||ai|| = ... ||a„|| = 1, the RHS can be written as 


ii6 - E ~ + E ^E + ^E j 

2=1 2=1 2=1 3^'i- ^ 

n n 

= \\b— 'y ( AiOiip — £(i + y ( Ai) 

i=l i=l 

= \\b-xf-£{i + \\x\\Ay 

Denote by a = ||a;|| and let p be the projection of b to conejai,..., a„}, which implies ||6 — p|| = 
min{||6 — a:|| | a: € conejai,..., a„}}. 


Claim. For all b, x, a,p given above, we have ||6 — a;|p > — 2a||p|| + 1. 


By this claim (proved later), we have: 


\\T{b)-T{x)f 
>c?- 2a||p|| + 1 - £(l + pAof 


> - 2a||p|| + 1 - e(l + II^IU)^ 

= (l-£MA)a^-2(|b||+eMA)a + (l-e)- 


The last expression can be viewed as a quadratic function with respect to a. We will prove this function 
is nonnegative for all a G K. This is equivalent to 


(Ibll + £dAf - (1 - eMa) (1 - e) < 0 

(ma + 2||p||ma + l)e < 1 — Ibll^ 

^ ,< 1-IHP ^_ 

-Mi + 2|bibA + i Mi + 2|HbA + i’ 


which holds for the choice of e as in the hypothesis. In summary, if the event 5'e occurs, then ||T(&) — 
T{x)\\‘^ > 0 for all x G conejai,..., a„}, i.e. T{x) ^ cone{T(ai),..., T(a„)}. Thus, 


Prob(r(6) i TC) > Prob(5'e) > 1 - 2n(n + 


as claimed. □ 

Proof of the claim. If x = 0 then the claim is trivially true, since ||6 — x||^ = ||6|p = 1 = — 2Q;|b|| + 1. 

Hence we assume x ^ 0. First consider the case p ^ 0. By Pythagoras’ theorem, we must have 
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cP = 1 — ||p|p. We denote by z = '^x, then ||z|| = ||p||. Set <5 = n^, we have 

\\b-xr = \\b-szr 

= {i-s)\\br+is^-s)\\zr+s\\b-zr 

= {i-6)+{6^-s)\\pr+s\\b-zr 

> {1 - 6) + {6^ - S)\\p\\^ + Sd^ 

= (l-<5) + (<5^-<5)|bf+ 5(l-|b|p) 

= 5^pr-26\\pr + l = a^-2a\\p\\ + l. 

Next, we consider the case p = 0. In this case we have b^{x) < 0 for all x G conejai,..., a„}. Indeed, 
for an arbitrary d > 0, 

0 < t(II^ “ bx\\^ — 1) = i(l + — 26b^x — 1) = — 2b"'"x 

0 0 

which tends to —2b^x when 5 —>■ O’*". Therefore 

||6-x||^ = 1 - 26^x + ||x||^ > ||x|p + I = - 2a||p|| + I, 

which proves the claim. □ □ 


6 Computational results 


Let T be the random projector, A the constraint matrix, b the RHS vector, and X either R" in the case of 
LP and Z" in the case of IP. We solve Ax = bAx G X and T{A)x = T{b)Ax G to compare accuracy and 
performance. In these results, A is dense. We generate (A, b) componentwise from three distributions: 
uniform on [0,1], exponential, gamma. For T, we only test the best-known type of projector matrix 
T(y) = Py, namely P is a random k xm matrix each component of which is independently drawn from a 
normal A/’(0, -^) distribution. All problems were solved using CPLEX 12.6 on an Intel i7 2.70GHz CPU 
with 16.0 GB RAM. All the computational experiments were carried out in JuMP [Sj. 

Because accuracy is guaranteed for feasible instances by Lemma 13.51 (i), we only test infeasible LP 
and IP feasibility instances. For every given size m x n of the constraint matrix, we generate 10 different 
instances, each of which is projected using 100 randomly generated projectors P. For each size, we 
compute the percentage of success, defined as an infeasible original problem being reduced to an infeasible 
projected problem. Performance is evaluated by recording the average user CPU time taken by CPLEX 
to solve the original problem, and, comparatively, the time taken to perform the matrix multiplication 
PA plus the time taken by CPLEX to solve the projected problem. 

In the above computational results, we only report the actual solver execution time (in the case 
of the original problem) and matrix multiplication plus solver execution (in the case of the projected 
problem). Lastly, although Tables [T] tell a successful story, we obtained less satisfactory results with 
other distributions. Sparse instances yielded accurate but poorly performant results. So far, this seems 
to be a good practical method for dense LP/IP. 
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Table 1: LP: above, IP: below. j4cc.: accuracy (% feas./infas. agreement), Orig.: original (CPU), Proj.: projected 
instances (CPU). 


m 

n 

Acc. 

Uniform 

Orig. 

Proj. 

Exponential 

Acc. Orig. Proj. 

Acc. 

Gamma 

Orig. 

Proj. 

600 

1000 

99.5% 

1.57s 

0.12s 

93.7% 

1.66s 

0.12s 

94.6% 

1.64s 

0.11s 

700 

1000 

99.5% 

2.39s 

0.12s 

92.8% 

2.19s 

0.12s 

93.1% 

2.15s 

0.11s 

800 

1000 

99.5% 

2.55s 

0.11s 

95.0% 

2.91s 

0.11s 

97.3% 

2.78 s 

0.11s 

900 

1000 

99.6% 

3.49s 

0.12s 

96.1% 

3.65s 

0.13s 

97.0% 

3.57s 

0.13s 

1000 

1500 

99.5% 

16.54s 

0.20s 

93.0% 

18.10s 

0.20s 

91.2% 

17.58s 

0.20s 

1200 

1500 

99.6% 

22.46s 

0.23s 

95.7% 

22.46s 

0.20s 

95.7% 

22.58s 

0.22s 

1400 

1500 

100% 

31.08s 

0.24s 

93.2% 

35.24s 

0.26s 

95.0% 

31.06s 

0.23s 

1500 

2000 

99.4% 

48.89s 

0.30s 

91.3% 

51.23s 

0.31s 

90.1% 

51.08s 

0.31 

1600 

2000 

99.8% 

58.36s 

0.34s 

93.8% 

58.87s 

0.34s 
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